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ABSTRACT 

We use an analytic model to study how inhomogeneous hydrogen reionization affects the temperature 
distribution of the intergalactic medium (IGM). During this process, the residual energy of each 
ionizing photon is deposited in the IGM as heat, increasing its temperature to 20, 000-30, 000 K; 
subsequent expansion of the Universe then cools the gas. Because reionization most likely proceeds 
from high to low densities, underdense voids are ionized last, have less time to cool, and are (on 
average) warmer than mean-density gas immediately after reionization is complete (an "inverted" 
equation of state). From this initial configuration, the low-density gas cools quickly and eventually 
returns to a more normal equation of state. The rapidly evolving temperature introduces systematic 
uncertainties in measurements of the ionizing background at z ~ 6. For example, late reionization 
implies rapid cooling, so that the ionizing background would have to evolve even more rapidly at 
z ~ 5-6 than typically claimed. This degeneracy is difficult to disentangle, because the Lya forest 
probes only a narrow range in densities (over which the gas is nearly isothermal). However, higher 
Lyman-series transitions probe wider density ranges, sampling different effective temperatures, and 
offer a new way to measure the IGM equation of state that should work where nearly saturated 
absorption precludes other methods. This will help to separate evolution in temperature from that 
in the ionizing background. While more detailed study with hydrodynamic simulations is needed, we 
show that such measurements could potentially distinguish early and late reionization using only a 
handful of lines of sight. 

Subject headings: cosmology: theory - intergalactic medium 



1. INTRODUCTION 

The two most dramatic events in the history of the 
intergalactic medium (IGM) are the reionization of hy- 
drogen (by the first generations of galaxies) and helium 
(by quasars). These have become key landmarks for 
both observational and theoretical cosmologists in the 
past several years. Evidence for hydrogen reionization 
comes from a number of directions, none of them clear 
but all consiste nt with (possibly extended) reionization 
at z ~ 6-10 (see lFan et al]|2006at IFurlanetto et aTll2006bl 
for recent reviews). Helium reionization is thought to 
occur during the quasar era, at z ~ 3, with a wide 
variety of supporting evidence - though much of it is 
controversial, and a clea r picture has yet to emerge (see 
iFurlanetto fc O h 2008b! for a recent summary of the ob- 
servations) . 

These two reionization epochs are largely responsi- 
ble for determining the thermal history of the IGM. 
Before hydrogen reionization, the neutral IGM cooled 
adiabatically until the first structures formed, proba- 
bly reaching temperatures T < 10 K. X-rays from 
the first galaxies most likely slowly heated the neutral 
IGM to T < 1000 K (fQhll2001b IVenkatesan et ai] |200H; 
iKuhlen et all 120061 : lFurlanettoH2006l ). However, hvdro- 
gen reionization caused a much more dramatic change: 
the ~ 5-10 eV leftover from eac h ionizing photon heated 
the IGM to ~ 2-3 x 10 4 K (iMiralda-Escude fc Reei 
1991 lAbel fc HaehneltJ Il999l : iTittlev fc Meiksinl 120071 : 
Trac et al.l 120081 )" The harder photons responsible for 
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helium reionization could have reheate d the IGM to 
simil ar or even larger , temp er atures (|Hui fc Gnedinl 
19971 IFurlanetto fc Obi l2008al iMcQuinn et alJ 120081 : 
Bolt on et alj|2008al) . " 

Once reionization is complete, this heating channel 
slows dramatically - because only the relatively small 
fraction of ions that recombine couple to the photoioniz- 
ing background. The subsequent temperature evolution 
is determined primarily by a balance of adiabatic heat- 
ing and/or cooling, photo-heating, and Compton coo l- 
ing (jMiralda-Escude fc Reeslll994j: iHui fc Gnedinlll997h . 
The competition between these processes forces the gas 
temperature to approach an asym ptotic form set by 
the background ioniz ing spectrum (|Hui fc Gnedinl [1997L 
IHui fc Hai man 2003). Because the magnitude (and in- 
deed sign) of the adiabatic term depends on whether the 
gas is over- or underdense, the IGM assumes an equa- 
tion of state T Tb(l + <5) 7_1 , where 8 is the fractional 
overdensity of the gas element and 7 > is nearly inde- 
pendent of 5 in most simple models. 

This equation of state inevitably affects many observ- 
ables of the Lya forest. In particular, the observed 
temperature evolution has been used to constrain the 
epochs of helium and hydrogen reionization: most dra- 
matically, the velocity widths of Lya for est absorbers 
seem to increase sharply at z ~ 3.2 (|Schave et alJ 

2000) , while the equation of state simult aneously flattens 
(jSchave et alJ 120001: IRicotti et al.ll2000h This may be a 
resul t of helium reionization (though seelMcDonal d etldl 

2001) . Recent models of inhomogeneous helium reion- 
ization show that these features are consistent with the 
behavior expected near that event dGleser et all 120051: 
IFurlanetto fc Obi l2008afc IMcQuinn et alJ l2008f ). More- 
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over, these models showed that helium reionization could 
dramatically transform the IGM thermal structure: if 
low-density voids are ionized last, they suffer the least 
amount of cooling and hence contain the hottest gas at 
the end of reionization, leading to an "inverted" equa- 
tion of state. Observations now suggest th at something 
simil a r may indeed be occ urring at z ~ 3 (|Becker et al.l 
120071 IBolton et al.l [2008bh . However, helium reioniza- 
tion is probably driven by rare, bright quasars; in that 
case, the effects are less dramatic because the nearly- 
random distribution of sources makes the process more 
stochas tic, and the equation of state may flatten but not 
invert (IFurlanetto fe Ohl l2008bt iMcQuinn et~ai1 120081 : 
IBolton et alj|2008a|) . 

Hydrogen reionization is an even more dramatic 
episode in the thermal history of the IGM, because it 
increases the IGM temperature by at least an order 
of magnitude. Moreover, it is most likely driven by 
large numbers of small galaxies, so the ionization topol- 
ogy (and hence thermal struct ure) are driven largely 
by the underlying den sity field ( Barkan a fc Loebl 12004 
iFurlanetto et al.ll2004f ). Unfortunately, helium reioniza- 
tion erases any information from this phase, so it is much 
more difficult to probe. One approach is to use the direct 
temperature measurements at z > 4 as a "fossil" record 
of hydrogen reionization. The inferred IGM temperature 
is relatively large, so this appears to require reionizatio n 
at z < 10 (|Theuns et all l2002d; iHui fc Haimanl 120031 ). 
However, because the IGM quickly approaches the afore- 
mentioned asymptotic equation of state, such inferences 
are difficult. 

Instead, it would be better to study the thermal effects 
of reionizatio n soon after that era ends, when the cooling 
is most rapid (|Trac et al.ll2008l ). This is, unfortunately, a 
difficult proposition, because the Lya forest - our princi- 
pal tool to study the IGM - becomes near l y satur ated in 
absorption at z > 5 (e.g.. iFan et~aT1l2001l . l2006bh . Thus 
no serious attempt has been made to measure the ther- 
mal structure at these redshifts. Nevertheless, it is cru- 
cial to interpreting the observations: the (temperature- 
dependent) recombination rate a(T) and the amplitude 
of the ionizing background T affect the (observed) op- 
tical depth only in the combination T/a(T), so inde- 
pendent estimates (or robust priors) on the tempera- 
ture are required to accurately constrain the growing 
ionizing background at z ~ 6 (|Fan et all 120021 [2006bt 
IBolton fc Haehnel"tll2007h . The rapid cooling expected 
after reionization ends may have serious implications for 
our understanding of this epoch. Moreover, recent sim- 
ulations show that the equation of state may be exceed- 
ingly complex shortly after reionization - with both an 
inver s ion and a large amount of scatter (|Bolton et al.l 
120041: iTrac et al.ll2008D . By contrast, almost all studies 
constraining T assume an isothermal IGM with no tem- 
perature evolution, or marginalize over simple equations 
of state. 

Indeed, if we could measure the evolution of the equa- 
tion of state as a function of redshift, we could place 
constraints on both the topology and redshift of reioniza- 
tion. The initial conditions (particularly the presence of 
an inverted equation of state) imprint information about 
the reionization topology, while the rate of evolution - 
which is most rapid immediately after reionization - tells 
us when that event occurred. Unfortunately, the satura- 



tion of the Lya forest, as well as the paucity of data at 
z ~ 6, imply that conventional techniques to measure the 
equation of state fail at high redshift. In this paper, we 
suggest that the flux transmission ratios of H I Lya and 
Ly/3, which probe different density regimes, can be used 
to measure the equation of state, provided the ionizing 
background does not vary systematically as a function of 
overdensity. This may open up studies of the equation 
of state at hitherto inaccessible redshifts, allowing us to 
probe its evolution over cosmic time. 

In this paper, we apply our model of inhomogeneous 
reionization (fully detailed in IFurlanetto fc OhT T2008a. 
hereafter FO08) to this important event. We will study 
how reionization affects the thermal structure of the 
IGM, and through that the Lya forest. We briefly de- 
scribe the model in We examine the resulting thermal 
histories in $3] and their effects on the Lya forest in 
We discuss the possibility of measuring the evolution of 
the equation of state in $51 We conclude in S|6] 

In our numerical calculations, we assume a cosmol- 
ogy with Q m = 0.26, Q A = 0.74, Q b = 0.044, H = 
100/ikms- 1 Mpc" 1 (with h = 0.74), n = 0.95, and 
as = 0-8, consistent w ith the most recent measurements 
(|Dunklev et al.ll200a) . Unless otherwise specified, we use 
comoving units for all distances. 

2. MODEL FOR THE THERMAL HISTORY 

Our model for inhomogeneous reionization is described 
fully in FO08. 3 In brief, it is based on the excursio n 
set approach to reionization of IFurlanetto et all |2004), 
which allows us to calculate the distribution of reion- 
ization redshifts for gas parcels as a function of their 
density. Given a final reionization redshift z r , the his- 
tory is determined by the assumed ionizing efficiencies of 
dark matter halos. Our fiducial model takes a constant 
ionizing efficiency in all halos with T v i r > 10 4 K (the 
threshold for atomic cooling in a hydrogen- helium gas), 
typi cal of many reionization models in the literature (see, 
e.g., iBarkana fc Loebll2001h . 

The excursion set model essentially evaluates the prob- 
ability that a small gas parcel of a specified density will 
lie inside an ionized bubble during various stages of reion- 
ization (and hence determines the local reionization red- 
shift for a given IGM patch). It implicitly assumes that 
the ionized bubbles trace dense regions surrounding the 
ionizing sources, so that dense gas elements are ionized 
relatively early, while voids must await the final stages 
of reionization. Note that, in FO08, we allowed for the 
possibility of "stochastic" reionization where the ioniza- 
tion history of a given gas parcel was independent of 
its density. While this is important for helium reion- 
ization, because that event is driven by rare and brigh t 
quasars (jFurlanetto fc Obi l2008bl : IMcQuinn et al.ll2008h . 
such a stochastic phase is probably not relevant during 
hydrogen reionization, which is driven by a multitude of 
faint galaxies. We will therefore focus exclusively on the 
density-driven model. 

Once the reionization history of an IGM parcel is cho- 
sen, we follow the su bsequent therm a l evol ution using 
a simplified version of IHui fc Gnedinl (|1997f ). which in- 
cludes all the relevant atomic heating, recombination, 

3 The original paper focused on helium reionization; here we 
only consider hydrogen. 
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Fig. 1. — Top: Median temperature at z = 6 and z = 5 (left and right panels, respectively) as a function of density after hydrogen 
reionization. The solid, long-dashed, short-dashed, dot-dashed, and dotted curves assume that reionization completes at z r = 6, 7, 8, 9, and 
10. Bottom: Local equation of state index, 7 — 1, for the same scenarios. All panels assume that the initial post-reionization temperature 
is T H = 30, 000 K. 



and photoionization processes as well as the density evo- 
lution (traced through the spherical collapse model; see 
FO08 for details). 4 After reionization, the dominant pro- 
cesses are photoionization heating and adiabatic cooling 
(with Compton cooling an important factor as well at 
the highest redshifts). 

2.1. The Initial Post- Reionization Temperature 

By construction, we cannot use our code while a gas 
parcel is first being ionized. It is during this brief pe- 
riod that the gas is rapidly heated by the leftover en- 
ergy of each ionizing photon. We instead set the initial 
post-reionization temperature, Th, by hand. A num- 
ber of factors that cannot self-consistently be included 
in our code affect this initial temperature - most im- 
portantly, we do not know the radiation spectrum that 
illuminates each gas element. This obviously depends on 
the UV spectra of the star-forming galaxies that provide 
the ionizing photons. The Starburst99 models show that 
low-metallicity galaxies have a luminosity density that 
varies roughly as L v cx v~ 2 bet ween the hydrogen an d 
helium ionization edges (Fig. 2e. iLeitherer et aLH l999). 
However, other factors are also involved - for example, 
spectral filtering by dense gas parcels (which may ab- 
sorb low-energy photons but allow high-energy photons 
to pass) makes the background spectru m inhomogeneous 
and h arder than the input spectrum (|Abel &: Haehneltl 
[19991) . 

We will therefore appeal to some general ar- 
guments that bound the likely value of Th 
(|Miralda-Escude fc Reed fl99i [AbeT fc Haehneltl fl999h . 
If the region between the gas element and the ionizing 

4 For computational simplicity, once a parcel is ionized we as- 
sume that it remains in ionization equilibrium with the radia- 
tion background. We have tested this assumption against a non- 
equilibrium code (see FO08) and find it to be an excellent approx- 
imation in all cases. 



sources is optically thin, the mean excess energy of ion- 
izing photons -Ethin is the average of the entire spectrum 
weighted by the ionization cross section, Oi cx T~ 3 , so 
low-metallicity stars yield (Ethin) /Ei ps 1/4, where Ei is 
the ionization potential. 5 However, the hydrogen density 
is so large that even most low-density gas systems will 
themselves be optically thick. In that limit, where all 
ionizing photons are absorbed, we do not weight by the 
cross section, so (£thick) /T"i ~ 3/5. In either case, this 
energy must then be shared with all the IGM baryons 
through Coulomb interactions; the net temperature 
change is then AT « 0.5(2/3/c B ) (E) ~ 30, 000 K for the 
optically thick case. In the following, we assume that 
Th = AT, which should be an excellent approxima- 
tion given the low IGM temperatures expected before 
reionization, eve n in the presence of X-ray heating (e.g., 
lFurlanettdl200l . 6 

Our code also requires that the ionizing background 
be specified once a gas parcel is ionized. For simplicity, 
we will take a constant power-law shape and amplitude 
at all redshifts, with the angle- averaged specific intensity 
J v cx v~ a from 13.6 to 54.4 eV, and a = 1.5 (correspond- 
ing to a slightly hardened stellar spectrum). In our cal- 
culations, we set the ionization rate T = 10~ 12 ri2 s _1 
to a relatively high value that is appropriate well after 
reionization, Ti2 = 1. Fortunately, these assumptions 
make almost no difference to our results. Although we 
do not allow the amplitude of the ionizing background 
to evolve (of course it actually increases as more and 
more of the IGM is ionized), this actually has no di- 

5 Here we have ignored absorption by neutral helium and in- 
cluded all photons with 13.6 eV < E < 54.4 eV; in reality, ~ 1/8 
of photons above 24.6 eV will ionize helium instead, producing less 
heat than our standard expression. But this provides only a small 
correction to our already-uncertain estimates. 

6 Note that FO08 assumed a significantly lower post-hydrogen 
reionization temperature, Th = 15, 000 K, and so found much 
weaker effects at 2 ~ 4 than we do here. 
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Fig. 2.— As Fig.[j] but for T H = 20,000 K. 



rect effect on the thermal history: a larger T increases 
the rate at which neutral atoms interact with the back- 
ground, but it reduces the neutral fraction by exactly the 
same amount. It only weakly affects secondary heating 
and cooling mechanisms, such as Compton scattering, by 
slightly increasing the density of free electrons. In addi- 
tion, note that, after the initial reionization, the shape 
of the ionizing background also has little effect on the 
temperatures. 

3. THE TEMPERATURE HISTORY 

The top panels of Figure [1] show the median IGM tem- 
perature as a function of density (here p/p — 1 + 5) for 
gas at z = 6 and 5 (left and right panels, respectively), 
if reionization occurs at z r — 6, 7, 8, 9, and 10 (from top 
to bottom in each panel) if T H = 30, 000 K. The bot- 
tom panels show the effective equation of state at each 
density, 



Figure^] shows similar information for models with Tjj = 
20,000 K. 

If reionization occurs early enough, the competition 
between photohcating and adiabatic cooling (including 
both the expansion of the universe and structure forma- 
tion, whose effect depends on 5) dr ives the IGM toward 
an asymptotic equat ion of state (|Hui k. Gnedinl 119971 ; 
iHui fc Haimanl I2003D ; at the high redshifts of interest 
here, this process is accelerated significantly by Comp- 
ton cooling. This asymptote is typically described as a 
simple power law, T oc (1 + (5) 7 " 1 , which we find to be 
a good approximation: reaching (7 — 1) ~ 0.6 long after 
reionization, with some flattening at very high densities. 
However, at redshifts closer to the time of reionization, a 
simple power law obviously no longer suffices to describe 
the results. Although 7 remains nearly the same at high 
densities, there is substantial flattening at lower densities 
- and, sufficiently close to reionization, it even becomes 
inverted (7 — 1 < 0) near and below the mean density. 
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This occurs because the highest-density regions - 
which preferentially lie near galaxies - are ionized first, 
while the lowest-density voids are far from sources and so 
get ionized last. All are initially heated to the same tem- 
perature (Th), but the former have longer to cool. (The 
highest density regions undergo rapid enough adiabatic 
heating that they stay warm.) Thus, when reionization 
ends, the hottest regions correspond to the most-recently 
ionized voids. Recent num erical simulations of reioniza- 
tion by iTrac et all (|2008l ) show a very similar median 
relation (although more similar to our Fig. [5] than the 
higher Tjj case). An inverted equation of state could 
also occur during helium reionization, if small sources are 
common enough (FO08), although in practice it appears 
unlikely given the dominance of b right quasars and the 
more complex radi ative transfer (|McQuinn et al.l l2008h 
iBolton et aLll2008aD . 

The IGM cools rapidly after this (especially with the 
addition of Compton cooling at high redshifts) , and most 
quickly in the voids, so this inversion does not last long. 
Once z < z r — 1, the low density gas has become nearly 
isothermal, from there slowly returning to a monotonic 
equation of state and eventually the thermal asymp- 
tote. Thus, photoheating will only substantially affect 
the equation of state for a relatively brief period after 
hydrogen reionization. In particular, Figure [3] shows that 
the effects will be relatively modest by z = 4: the differ- 
ences at the mean density are only ~ 2000 K, although 
they remain a factor of a few at the lowest densities. 

There have been attempts to c onstrain z r from the 
Lyq forest temperatu re at z ~ 4 (|Theuns et alj [2002c: 
IHui fc Haimanl 120031 ). under the assumption that the 
temperature is uncontaminated by helium reionization 
there. Unfortunately, the data are not nearly good 
enough to measure the relatively small expected dif- 
ference. In fact, temperatures inferred from the exist- 
ing data generally lie well above predictions that in- 
clude only hydrogen p hotoheating (|Schave et al.1 120001 : 
iZaldarriaea et al.ll2001f ). 
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Figure H shows the second major effect of inhomoge- 
neous reionization on the IGM temperature distribution: 
because different gas parcels (even at the same local den- 
sity) have different ionization histories, reionization cre- 
ates a significant scatter in the temperature distribution. 
The curves in each panel show the 10th, 25th, 50th, 75th, 
and 90th percentiles of the distributions, from bottom to 
top; the different panels show different times after reion- 
ization, all assuming z r = 6. 

Again, because reionization proceeds "inside-out," the 
highest density regions are ionized first, with relatively 
little scatter in their reionization rcdshifts (which is then 
decreased even further by adiabatic heating). The low- 
est density regions are ionized last, again with relatively 
little scatter: in the excursion set formalism, such voids 
are always near the last to be ionized - in this case at 
6 < z < 6.5 - so the final spread in temperatures is mod- 
est. However, the probability that gas near the mean 
density has been ionized at any time in the central pe- 
riod of reionization is roughly equal to the mean ion- 
ized fraction, so there is a considerable spread in the 
temperatures around S ~ 0. The factor of two spread 
persists until z ~ 5, but the relation then tightens as 
gas approaches the thermal asymptote, and by z = 4 
the scatter is small. In principle, this reduction in scat- 
ter could be a tell-tale signature of when reionization 
took place, as has bee n suggested for helium reionization 
(|Hui fc Haimanl 12003ft . although that has not yet bee n 
detected (e.g., iTheuns et all l2002bt IZaldarriagal l2002f K 
For hydrogen reionization, it is unfortunately difficult to 
obtain sufficient statistics to measure the scatter in the 
equation of state at high redshifts, while the lower red- 
shift observations at z ~ 4 are likely contaminated by 
helium reionization. 

The reionization simulations of lTrac et al.l (|2008l ) show 
comparable scatter near the mean density but much more 
in low-density voids and along the high-density tail. This 
suggests that our excursion set model relating the large- 
scale ionized regions to small gas parcels may contain 



too much of a correlation across scales, although there 
are other complexities in the simulation (in particular, 
differing source prescriptions, radiative transfer, and the 
complex geometry of the cosmic web) that can explain 
these differences. In any case, increased scatter is likely 
to strengthen the observable consequences that we de- 
scribe below. 

4. EFFECTS ON THE HI Lya FOREST 

We have seen in the previous section that the 
equation of state evolves significantly after reioniza- 
tion completes. Here we will consider the impli- 
cations of this phenomenon for observables, and in 
particular the HI Lya forest. During the helium 
reionization era, the corresponding effects are sub- 
stantial: the temperature increase and flattening of 
the equation of state may have already been ob- 
served directly, through the line width distribution in 
the L ya forest dSchave et all I1999L l2000t iRicotti et all 
120001 though see also iMcDonald et al.l l2001h . and in 
directly, p ossibly through the e volving effective opti- 
cal d epth (|Bernardi et al.l 120031 : iFaucher-Giguere et al.l 
2008). The observed behavior matches the qualitative 
expectations of theoretical studies of helium reioniza- 
tion, although not n ecessarily the details - especially 
for th e optical depth ([Theuns et alj|2002at iGleser et al.l 
120051: iFurlanetto fc Obi l2008al iMcQuinn et alJ 120081 : 
iBolton et al1l2009f K 

However, extracting the details of the equation of 
state following hydrogen reionization promises to be 
much more difficult. Direct methods, which rely on the 
linewidths of the features, are useless because the for- 
est is so saturated - with individual features no longer 
identifiable. Instead, we will likely need to constrain the 
thermal evolution indirectly, through the overall trans- 
mission in the forest. 

To examine how the optical depth evolves, we fol- 
low the model for the I GM density structure of 
iMiralda-Escude et alJ (|2000T ). They found that a rea- 
sonable fit to the volume-averaged density distribution 
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Fig. 4. — Distribution of IGM temperature as a function of density, during and after rcionization. All panels assume z r = 6 and take 
z = 4, 5, 5.5, and 6, clockwise from bottom left. The curves show the 10th, 25th, 50th, 75th, and 90th percentiles of the distributions. 
Both panels assume that the post-reionization temperature is Tn = 30, 000 K. 



of IGM gas, fV(A) (where A 
z ~ 2-4 is: 

Py(A) dA = AoA^exp 



1 + <5), in simulations at 



(A" 



-2/3 



Co) 



2(2 W (2) 
Intuitively, the underlying Gaussian density fluctuations 
are modified through nonlinear void growth and a power 
law tail at large A. Here <5o essentially represents the 
variance of density fluctuations smoothed on the Jeans 
scale for an ionized medium; thus 5q cx (1 + z)^ 1 at 
high redshifts (such as those here). The power-law ex- 
ponent (3 determines the behavior at large densities; for 
isothermal spheres, it is ft = 2.5, which we assume to be 
appropriate at z = 6. One c aveat is necessary here: the 
Miralda-Escude et al. (2000) simulations included Jeans 
smoothing from a specific reionization history, so we can- 
not use it to describe the dynamical evo lution of the 
gas d istribution that follows reionization (|Pawlik et al.l 
12001 . 

We will only consider the most easily measured char- 
acteristic of the absorption: the effective optical depth, 
r e ff, which simply parameterizes the mean transmission 
in Lya, T a = exp(— T e ff ia ). We estimate the transmission 
by simply integrating over T and A, 



dAPv-(A) / dTp(T\A)e 



-t«(A,T) 



(3) 



where p{T\ A) is the probability that a gas parcel of den- 
sity A has a temperature T and is given by our inhomoge- 
neous reionization model. With the " fluctuating Gunn- 
Peter son trough approximation" (e.g-.lGunn fc Peterson! 
119651: iRauch et al.lll997h ICroft et al.lll998f) . the Lyq op- 
tical depth of a gas parcel is approximately 

A 2 T -0.7 / 1 + z \ 9 / 2 

Tq (A,T)«13^- (1±*\ , (4) 
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where we have used the case- A recombination coefficient, 
ola ~ 4.2 x 10 _13 T 4 ~ 7 cm 3 s^ 1 (appropriate for the low- 
density diffuse IGM, where most of the transmission will 



occur), T4 — T/(10 4 K), and we have assumed a large 
enough redshift that the cosmological constant can be ig- 
nored. Unfortunately, this procedure does not exactly re- 
produce the features of the Lya forest, and it is necessary 
to renormalize r a in order to match measure ments from 
more detailed studies of the forest. Following iFan et al.l 
(2002) we introduce a correction factor k — 0.3 on the 
right hand side of equation Q t o reproduce the T e s. a —T 
relatio n at z = 4 found by [McDonald & Mira lda-Escudel 
(2001). Although this is a large factor, the hope is that 
it does not evolve significantly with redshift, so that our 
results at other redshifts are reasonable. 

The crucial point is that r Q tx a(T)/r, so the am- 
plitude of the ionizing background cannot be measured 
independently of the temperature. In practice, given 
the difficulty of independently measuring T at the high 
redshifts of interest for reionization, studies typically 
take it to be co n stant b oth with time and density (e.g., 
IFan et all 120011 . [2006b) or margina lize over it with a 
range of simple equations of state ([Bolton &: Haehneitl 
2007). Our models show that this is clearly too simple: 
the Lya forest is most sensitive to A ~ 0.2 at z ~ 6, 
where the temperature can vary by nearly an order of 
magnitude after reionization (see Fig. [J). Naively, this 
leads to uncertainties in T of a factor of a few. Even 
more worrying are the systematic errors introduced by 
post-reionization cooling, which affect estimates of the 
evolution of L. 

The filled symbols in Figure Ogive a practical example 
of this thermal evolution's effects on T c ff iQ . Here we show 
how the transmission evolves for three reionization sce- 
narios (with z r = 6, 7, and 10 for the triangles, squares, 
and hexagons). For concreteness, we take Ti2 = 0.1, 0.2, 
and 0.4 at z — 6, 5, and 4 (chosen to match roughly 
the estimates from lFan et al.l [20 06b). The red and blue 
symbols take T H = 20, 000 and 30, 000 K, respectively. 

As expected, the differences between the three scenar- 
ios are small at z = 4, because all the models give nearly 
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Fig. 5. — Top panel: Effective optical depths for Lya (filled 
symbols) and Ly/3 (open symbols) in three different reionization 
scenarios: z r = 6, 7, and 10 (triangles, squares, and hexagons, 
respectively). In all cases, Ti2 = 0.1, 0.2, and 0.4 at z = 6, 5, and 
4, respectively. The red and blue symbols take Th = 20, 000 and 
30, 000 K, respectively (note that they are slightly displaced from 
2 = 6, 5, and 4 only for plotting purposes). Bottom panel: Ratio 
of effective temperatures T a /Tp = (A a /Ap) 0,7 , for the models 
with z r = 7 and 10 (solid and dotted curves, respectively), with 
Tjj = 30, 000 K, depicted in the top panel. 



the same temperature at that point (the difference is ~ 5 
per cent in absolute transmission). However, as we ap- 
proach z = 6, the transmission increases dramatically 
for the z r = 6 model: it is over an order of magnitude 
larger than that of the z r — 10 model, purely because 
of the high temperatures achieved in low-density voids 
just after reionization. This higher temperature helps to 
compensate for the increasing density of the Universe; 
in other words, the rapid cooling following reionization 
slows the evolution of T e s,a- 

The open symbols in Figure O show the same quan- 
tity, but for the Ly/3 transition. This line has a signifi- 
cantly smaller oscillator strength (leading to a prefactor 
6.24 times smaller in eq. H}; it therefore samples higher 
densities, whose temperatures remain more stable after 
reionization. As a result, r c g^ varies by a much smaller 
amount across these reionization models (differing by a 
factor of only ~ 2.5 in absolute transmission at z = 6). 

Figure [6] provides some intuition about this behavior. 
The left panel shows the logarithmic contributions by 
density to the outer integral in equation ([3]), or in other 
words the range of densities that are actually visible in 
the forest. The black curves are for Lya while the blue 
curves are for Ly/3. Within each set, the solid, dashed, 
and dotted curves assume that reionization completes at 
z r = 6, 7, and 10, respectively. All take T H = 30, 000 K 
and the same ionizing backgrounds as in Fig. [5] The 
right panel shows the median overdensity contributing to 
the total IGM transmission, as well as the overdensitics 
between which 68% of the flux transmission occurs. The 
figure assumes an isother mal equa t ion of state and uses 
a spline fit to the data of ISongaiTal ([2004) for T a . 

Several key points are obvious from the figure. First, in 



the left panel all the curves converge at A < 0.15. This is 
becau se voids with A < 0.2 a r e exp onentially suppressed 
in the iMiralda-Escude et alj (|2000h density distribution 
and is a phenomenon independent of our temperature 
distributions. On the other hand, the turnover at high 
densities is primarily due to the e~ T factor and depends 
on both the transition and the temperature. For exam- 
ple, the dot-dashed curve in (a) shows APy(A) at z = 4. 
All the other curves trace it closely until r ~ 1. A second 
feature is also crucial: at z = 6 the largest transmission 
occurs when z r — 6. The order-of-magnitude increase 
in temperature opens up a substantial range of densi- 
ties (over which Py is increasing rapidly), raising the 
transmission by a large factor. These points are ampli- 
fied by the right panel, which shows the density depen- 
dence with redshift, given the the empirically observed 
Lya flux transmission. While the median overdensity for 
Lya and Ly/3 transmission do not differ strongly, the lat- 
ter transition accesses a wider range of densities above 
the median. On the other hand, the range of contri- 
butions for both transitions narrows significantly toward 
high redshift (when increasing densities and a weaker ion- 
izing background imply that higher overdensities do not 
transmit flux - leaving the observations heavily weighted 
toward voids). 

5. CAN WE MEASURE THE IGM TEMPERATURE AT 

Z > 5? 

From a practical standpoint, the most important goal 
is to break the degeneracy between T and the thermal 
structure. In one sense, the narrow range of densities 
probed by the Lya forest is helpful, because it insulates 
us from uncertainties in the equation of state parameter 
7: the transmission in all of our models can be closely 
mimicked by choosing an isothermal equation of state 
with appropriate normalization T. On the other hand, 
this same narrow range makes it even more difficult to 
measure evolution in the temperature, which shows up 
most clearly in the steepening of the equation of state 
with time. 

However, there is some hope, thanks to the higher- 
order Lyman transitions. Figure [S] shows that Ly/3 is 
sensitive to a wider range of densities than Lya, so it 
will effectively have a higher temperature at any redshift 
(and Ly7, with an optical depth 18 times smaller than 
Lya, will sample a still larger range). Intuitively, measur- 
ing the effective temperature at several different densities 
through these different transitions provides a handle on 
the shape of the equation of state, which can then be 
combined with models like our own to break the T-T de- 
generacy (and even to measure the reioniz ation redshift 
from the evolving thermal p roperties, as in lTheuns et alJ 
[2002cllHui & Haimarij200aL 

iDiikstra et al.l ( 2004 ) first recognized that the Ly/3 for- 
est could be used to measure the equation of state, be- 
cause it probes the IGM at higher overdensities rela- 
tive to Lya absorption. In particular, they showed that 
the small-scale Ly/3 power spectrum was sensitive to the 
equation of state, whereas the Lya power spectrum was 
relatively robust to such variations, due to its narrow 
density coverage. However, for our purposes it is not 
possible to measure the Lya or Ly/3 power spectra at 
high-redshift, given the paucity of data and the frequency 
of saturated absorption. Instead, we point out that the 




Fig. 6. — Left panel: Contribution by density to the total IGM transmission at z = 4 (a) and z = 6 (6). The black curves are for Lyo 
while the blue curves are for Ly/3. Within each set, the solid, dashed, and dotted curves assume reionization completes at z r = 6, 7, and 10, 
respectively. All curves assume Th = 30,000 K and T\2 = 0.1 [z = 6) and 0.4 (z = 4). In the top panel, the dot-dashed curve shows the 
density distribution from eq. (O- Right panel: Solid lines show the median overdensity which contributes to the total IGM transmission 
for the respective transitions (Lya and Ly/3, bottom and top), while the dashed (dotted) lines show the overdensities between which 68% 
of the flux transmission occurs for the respective transitions. 



relative effective optical depths in Lya and Ly/3 absorp- 
tion can already give us information about the equation 
of state. 

More formally, we have T e &,i = T e s } i(Ai) for each tran- 
sition i = a, f3, and 7, where A, = T/a(Ti) and Ti is 
the effective temperature for that transition. As de- 
scribed above, these Tj will differ because each transi- 
tion probes a different range of densities. On the other 
ha nd, T is (naively) indepe ndent of the transition. As 
in I Fan et all (|2001l I2006H ). we can then extract the 
Ai and measure the slope of the equation of state via 
A a /A p = (Ta/Tjg) - 7 . The redshift evolution of the slope 
can then reveal details of reionization and (after further 
modeling) T. For example, the bottom panel of Fig- 
ure [S] shows the inferred ratio in two reionization scenar- 
ios, with z r — 7 and z r = 10. Long after reionization, 
the temperatures no longer evolve, so T a /Tp is nearly 
constant (as in the z r = 10 scenario). But just after 
reionization, when the voids are cooling quickly, the ra- 
tio undergoes substantial evolution. It is this effect that 
can constrain the timing of reionization. 

For example, consider the line of sight to J1148+5251 
(a z = 6.42 quasar), which sho ws transmiss i on in 
the Lya, Ly/3, and Ly^ troughs (IWhite et all 120031: 
lOh fc Furlanettol [20051 : IWhite et all 120051 ). iFan et alJ 
(20063) used a technique very similar to ours, but as- 
suming an isothermal equation of state with T — 10 4 K, 
to transform the amount of transmission into measure- 
ments of r 12 = 0.06, 0.04, and 0.015 from the Lya, Ly/3, 
and Ly7 troughs, respectively. Although the accompa- 
nying errors are large (and difficult to estimate) the dis- 
crepancies are certainly worrisome. 

We have argued that these should instead be consid- 
ered measurements of FT®- 7 , so that (assuming a con- 
stant r) the three estimates actually imply T a w 2Tp w 
7T 7 - a strongly inverted equation of state; indeed, this 



requires a much steeper density dependence (7— 1 ~ —2) 
than seems reasonable (c.f. Fig. [I}. Preliminary calcula- 
tions for other quasars in the z ~ 4— 6 range also indicate 
that the majority of them have strongly inverted appar- 
ent equations of state. Clearly - at least according to our 
models - the post-reionization equation of state can ac- 
count for only a small fraction of these differences. Other 
effects may then be at play. 

For example, T may be much larger in low-density re- 
gions than in those probed by the higher-order transi- 
tions. At first blush, this seems surprising, as the source 
density inside voids is much smaller than that in dense 
regions, so that r Q /r^ < 1, making our derived T a /Tp 
lower limits. Of course, radiative transfer effects may 
help. But the reduced abundance of Lyman limit ab- 
sorbers in voids (which could increase the mean free 
path of ionizing photons and hence increase the radiation 
field) in unlikely to overcome the reduced abundance of 
galaxies (which are presumably more highly biased and 
hence more sensitive to the large scale density), in or- 
der to allow a net increase in T. So one must appeal 
to more complex effects, such as shadowing. For exam- 
ple, in the contex t of th e helium-ionizing background, 
iMaselli fc Ferraral (|2005h found that radiative transfer 
hardens the background in low-density regions. Whether 
a systematic increase in the ionizing background at lower 
overdensities can possibly account for the discrepancies 
remains unclear for now. 

A comparison to the simulations of iTrac et~a l. (2008) 
suggests another possibility: increased scatter in the 
equation of state (particularly if it is density-dependent). 
For a given mean equation of state T(A), scatter in- 
creases flux transmission, much as an inhomogeneous 
density distribution does. However, the increase is 
greater for T a than for Tp, implying that T e ff jQ /T e ff )( 3 will 
be biased downward, and hence (T a /Tp) will be biased 
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upward. This is because scatter in the equation of state 
causes flux transmission to be less heavily weighted to- 
ward rare voids (since hot denser regions become just as 
transparent). The steepness of Py(A) in the Lya trans- 
mission range (see Fig. |6|) then implies a large increase 
in the total transmission. Fluctuations in the ionizing 
background would have a similar effect. 

For instance, consider an isothermal mean equation 
of state (T and T are independent of density), where 
A = r_i 2 (_T/20,000K) - 7 has a lognormal distribu- 
tion with (A,ai n A) — (0.02,1), either due to density- 
independent scatter in the temperature or in the ioniz- 
ing background. This gives T e s j0l /T e s,f3 = 1-9 (as opposed 
to T e s t a/T c ff t i3 — 2.6 for an isothermal equation of state 
without scatter), i.e. increased flux transmission in Lya 
relative to Ly/3. We would then infer (T a /Tp) = 1.4 from 
our naive approach - an "inverted" equation of state - 
even though the gas is actually isothermal (in the mean) . 
Such biases will be further amplifie d if the sca t ter in - 
creases toward lower densities, as in iTrac et all (|2008| ). 
Note, however, that a full evaluation of the impact of 
scatter requires numerical simulations, because some of 
the scatter in the low density gas may not be important. 
For example, the high temperatures in low-density gas 
near large scale overdensities may not manifest them- 
selves in observations, since the absorption in such re- 
gions of the spectra is likely to be saturated anyway. 

Interestingly, both scatter in the equation of state 
and any density dependence in the ionizing background 
should fall with time (the latter because the photon mean 
free path increases with time), and so all of these effects 
only increase the likelihood of an "inverted" equation of 
state immediately after reionization, followed by a recov- 
ery toward the asymptotic relation. 

The above uncertainties notwithstanding, high-quality 
spectra of the sort already available might still provide 
some useful constraints. As an example, suppose that 
we have measured r c g^ a = 6.1 at z = 6. Models with 
z r = 6 and T 12 = 0.0385 or z r = 10 and r 12 = 0.1 both 
provide this level of Lya absorption, but 7"eflf„a = 2.34 
and 2.07 for the respective scenarios (or a ~ 20% differ- 
ence in transmitted flux). How many lines of sight are 
required to distinguish between these predictions? There 
are four sources of uncertainty on T e g t p: cosmic variance 
(which provides a maximal fractional error ~ 0.4 on the 
Ly/3 transmission for this value of opacity; iLidz et all 
I2006I ). 7 uncertaint y in the foreground L ya absorption 
( Tfg = 2.38 ± 0.32: foh fc Furlanettdl2005f ). uncertainties 
in estimating the quas ar continuum Fq (~ 10% at Ly/3; 
lOh fc Furlanettol 120 05) . and measurement error (which 
is much smaller tha n the other uncert ainties for high- 
quality spectra; e.g.. lWhite et al.ll2003f ). Thus the total 
uncertainty in the measured flux -Fobs over TV lines of 
sight is 



5F oha *^J(6T fs )* + (6T^+(^y (5) 

7 Actually, this value corresponds to the cosmic variance on dif- 
ferent lines of sight. We are comparing the Lya and Ly/3 forests 
along the same line of sight, so the cosmic variance should be much 
smaller - sourced only by the decorrelation between Lya and Ly/3 
from sampling different densities. This aspect must be calibrated 
with simulations. 



where the second term in the square root conservatively 
accounts for cosmic variance. Thus, only ~ 5 (20) lines of 
sight would be needed to distinguish these two reioniza- 
tion scenarios at 68% (95%) confidence - provided that 
the Lya transmission can itself be cleanly estimated (and 
provided we can robustly model the spatial distribution 
of r and the temperature scatter). Working with Lya is 
somewhat harder as absorption is more easily saturated, 
yielding only lower bounds on T c ff iCC , and hence upper 
bounds on T a /Tp. However, our results show that, with 
improved modeling, constraints on the timing of reion- 
ization should be possible. 

6. DISCUSSION 

We have described how inhomogeneous hydrogen 
reionization affects the thermal structure of the IGM. Be- 
cause it proceeds "inside out," f rom high to low densit y 
regions ([Barkana fc Loebl l2004t iFurlanetto et all l2004h , 
voids are ionized last and so are hot ter than denser gas 
near the end of reionization (see also lBolton et al.ll2004l : 
ITrac et al"1l2008l ). Rapid adiabatic and Compton cooling 
work quickly to erase this inversion, and the IGM equa- 
tion of state passes through a nearly isothermal phase 
(at A < 1) and then approaches its relatively steep and 
monotonic asymptote. Late reionization at z r = 6 leaves 
the low-density gas in the isothermal phase at z ~ 4; if 
reionization is earlier, it is close to the asymptotic form. 
In our model, the scatter around the median tempera- 
ture at low and high densities is relatively small, because 
extreme densities are all either ionized very near the end 
of reionization (for A <C 1) or very near the beginning. 
The scatter is larger near the mean density, because the 
spread in reionizat ion redshifts is m uch larger. The re- 
cent simulations of ITrac et al.l (|2008l ) show larger scatter 
at high and low densities, perhaps indicating that radia- 
tive transfer plays a significant role in the temperature 
distribution. 

Although similar trends may occur during helium 
reionization at z ~ 3 (FO08), "inside out" reionization 
is less relevant in practice during that phase, because 
the ionizing sources are so rare and bright that the cor- 
relation between large-scale overdensities and ionized re- 
gions is much weaker and because ra diative transfer is 
much more complex during this era (IFurlanetto fc Obi 
I2008bl iMcQ.uinn et alj|2008t iBolton et alll2008al) . From 
a theoretical standpoint, an "inside out" picture is much 
more likely to be accurate during hydrogen reionization, 
which is driven by large groupin gs of small sources whos e 
random fluctuations are small (jFurlanetto et al.ll2006al) . 

The Lya forest is so thick at z > 5 that individual lines 
can no longer be distinguished, so the techniques used to 
measure the IGM temperature directly at lower redshifts 
dSchave et al.ll2000t [Ricotti et all [200(1 iMcDonald et all 
I200I " are no longer useful. However, the thermal struc- 
ture does still affect the forest: we found that shifting the 
end of reionization over the range z r = 6-10 can change 
the total transmission by up to an order of magnitude at 
z = 6, although by z — 4 these differences have nearly 
disappeared (thanks largely to the decreasing mean den- 
sity of the IGM, which makes a much larger fraction of 
the IGM visible). 

This raises two important questions for understanding 
the high- z IGM. First, can the temperature evolution be 
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measured and used to constrain the timing of reioniza- 
tion? We have shown that the expected variation in tem- 
perature is modest at z < 4, so using that era is difficult 
(even leavin g aside possible contamination from helium 
reionization: iTheuns et aT1l2002d : iHui fc Haimanll2003f ). 
Moreover, observations yield significantly higher temper- 
atures than any published model, with error s much larger 
than the differences between the models (jSchave et al.l 
120001: iZaldarriaea et all 120011 ). At higher redshifts, de- 
generacy with the amplitude of the ionizing background 
is troublesome: an isothermal IGM is completely de- 
generate with r, and, for the Lya transition at z > 4, 
isothermality is typically an excellent approximation, be- 
cause the forest is only able to sample a very limited 
range in densities. 

However, higher-order transitions can help signifi- 
cantly (at least in principle), because they sample a wider 
range of densities (see Fig. [6]) and hence a different por- 
tion of the equation of state. Deviations from isother- 
mality can be measured via the transmission ratios and 
used to map the equation of state (and hence reionization 
rcdshift). This is complicated by the rapidly diminishing 
transmission as well as contamination from the lower- 
redshift Lya forest, but we have found that differentiat- 
ing z r — 6 and z r — 10 via the Ly/3 transmission only 
requires a modest increase in the number of known lines 
of sight to distant quasars. However, we will certainly 
need more sophisticated models that include variations 
in the ionizing background across these different environ- 
ments, as well as scatter in the equation of state. 

The second consequence is for measurements of T, 
which are particularly interesting near reionization. To 
date, claimed constraints on r at z — 6 have ignored 
thermal evolution, eith er assuming isothermality (e.g., 
iFan et al.l [20011 l2006bD or marginalizin g over standard 
equations of state (jBolton fc Haehneltl 12007). Here we 
have shown that the evolution is much more complex, 
which will systematically affect measurements of the ion- 
izing background. Unfortunately, with existing data it is 
not at all clear how these effects will play out. 

As shown in Figure [51 just after reionization both the 
Lya and Ly/3 optical depths evolve quite slowly, even if 
r is increasing, because rapid IGM cooling cancels out 
much of the change in T as well as the overall density 
evolution. If reionization occurs late, T must therefore 
evolve even more quickly than commonly assumed over 
Z ~ 5-6: in our models, if reionization ends at z = 6, so 
that (T a ) ~ 25, 000 K, we would require Ti 2 « 0.04 to at- 
tain T e s }a w 6 (as opposed to Ti 2 = 0.1 if z r = 10). Thus, 
if the conventional explanation is correct, the ionizing 
background would have to increase by a factor of five 
from z = 6 to z = 5; whether this is reasonable depends 
entirely on the presumed evoluti on of Lyman-limit sys- 
tems , which is highly uncertain (jFurlanetto fc Mesingerl 
2008). If instead reionization - or at least reheating - 
ends much earlier, T could increase relatively slowly over 
the same time interval. 

Without more detailed modeling, for now it is per- 
haps better to consider uncertainties in the equation 
of state to be "noise" in measurements of the ioniz- 
ing background. Our model suggests that such system- 
atic uncertainties can be relatively large, even compared 
to the substantial statis tical errors on the ex isting ob- 
servations; for example, iBolton et al.l (j2008al ) estimate 



ri 2 (T Q /10,000 K) 7 = 0.2±°j, 5 at z = 6. Allowing (T a ) 
to vary from 4, 000-25, 000 K increases the allowed range 
to Ti2 = 0-0.6. Thus, even today's constraints need to 
account for the thermal effects of reionization. 

It is important to note that our model neglects (at 
least) three aspects of photoheating that will affect the 
Lya forest, so more detailed simulations will be required 
for careful comparison to observations. First, our fluc- 
tuating Gunn-Peterson approximation (used in $4]) ig~ 
nores thermal broadening and does not perfectly repro- 
duce the line structure of the forest (jBolton et alj 120051 
iTvtler et all 12004 I Jena et al J [20051 ). For example, the 
heating that accompanies reionization will make the for- 
est lines broader, decreasing the number of saturated 
lines and increasing the overall optical depth. This will 
oppose the direct effect of heating (which reduces the 
neutral fraction) and so partially compensate for the de- 
creased optical depth following reionization. As the Uni- 
verse cools, the lines will become narrower and more sat- 
urated, again opposing the slow decrease in the overall 
transmission in our results. This will make photoheating 
even more difficult to observe. 

Second, we have ignored radiative transfer effects 
during reionization, where the radiation field hard- 
ens as it propagates o utward from an ionizing source 
(|Abel fc Haehneltl IT999) . Note that such spectral hard- 
ening effects could also lead to a systematic tilt in the 
initial equation of state, in that voids may see a more 
heavily filtered and hardened radia tion field, and henc e 
be heated to higher temperatures (jBolton et al.l [2004). 
Such effects may stymie efforts to infer the topology of 
reionization from an inverted equation of state. However, 
we expect them to be much weaker for hydrogen reion- 
ization as compared to helium reionization, since the for- 
mer case is driven by a softer (primarily stellar) radiation 
field. Furthermore, the shape of the UV radiation field 
does not significantly affect the post-reionization evolu- 
tion of the equation of state: the relatively low recombi- 
nation rate of hydrogen means that the gas largely loses 
thermal contact with the radiation field. Thus, efforts to 
infer the rcdshift of reionization using the techniques we 
have discussed may still be robust to radiative transfer 
effects, although this requires more detailed study with 
numerical simulations. Other tha n the increased sc atter, 
comparison to the simulations of iTrac et alj ([2008) docs 
not reveal any substantial biases in our approach. 

The final effect is more subtle: the increased temper- 
ature following reionization will also increase the pres- 
sure of the IGM, leading to a rearrangement of IGM 
gas through Jeans smoothing. While well-understood 
in the contex t of minihalos (where it is called "p hoto- 
evaporation" : IShapiro etal] 12004 Uliev et al J 120051) . the 
analog ous IGM processes have been considered only re- 
cently (jPawlik et al.ll2008T ). Jeans smoothing will move 
material away from high-density regions, "puffing out" 
filaments to fill voids and making Py(A) more s t rongly 
peaked around the mean density. iPawlik et all (2008) 
found that this adjustment occurs over approximately 
an expansion time (roughly the sound-crossing time at 
densities near the cosmic mean). This will affect both 
the mean transmission [through Py(A)] and the mean 
free pat h of ionizing photons (w hich depends on dense 
clumps; iFurlanetto fc Of] [2005) . However, our simple 
analytic model does not suffice to describe this process, 
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so we defer its inclusion to future, more detailed work. 

Nonetheless, it is worth noting that photoheating re- 
duces transmission at late times once these hydrody- 
namic effects set in, making the voids, which are the pri- 
mary source of transmission, higher density and so more 
opaque. As redshift decreases, this creates a systematic 
bias toward lower and lower estimates of T compared to 
the true value, just like the temperature evolution ef- 
fects we have been discussing. If so, T must evolve even 
faster than we have quoted above. Such effects also re- 
duce T a /Tp (since T a is more heavily weighted toward 
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